System size resonance in coupled noisy systems and in the Ising model 
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We consider an ensemble of coupled nonlinear noisy oscillators demonstrating in the thermody- 
namic limit an Ising-type transition. In the ordered phase and for finite ensembles stochastic flips 
of the mean field are observed with the rate depending on the ensemble size. When a small periodic 
force acts on the ensemble, the linear response of the system has a maximum at a certain system 
size, similar to the stochastic resonance phenomenon. We demonstrate this effect of system size 
resonance for different types of noisy oscillators and for different ensembles - lattices with nearest 
neighbors coupling and globally coupled populations. The Ising model is also shown to demonstrate 
the system size resonance. 

PACS numbers: 05.40.Ca, 05.45.-a, 05.50+q 
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Stochastic resonance has attracted much interest re- 
cently 0. As was demonstrated in a response of a 
noisy nonlinear system to a periodic forcing can exhibit a 
resonance-like dependence on the noise intensity. In other 
words, there exists a "resonant" noise intensity at which 
the response to a periodic force is maximally ordered. 
Stochastic resonance has been observed in numerous ex- 
periments Noteworthy, the order in a noise-driven 
system can have a maximum at a certain noise level even 
in the absence of periodic forcing, this phenomenon being 
called coherence resonance Q. 

Being first discussed in the context of a simple bistable 
model, stochastic resonance has been also studied in 
complex systems consisting of many elementary bistable 
cells f§, moreover, the resonance may be enhanced due 
to coupling Q . In this paper we discuss another type of 
resonance in such systems, namely the system size res- 
onance, when the dynamics is maximally ordered at a 
certain number of interacting subsystems. Contrary to 
previous reports of array-enhanced stochastic resonance 
(cf. also |7)), here we fix the noise strength, coupling, 
and other parameters; only the the size of the ensemble 
changes. 

The basic model to be considered below is the ensemble 
of noise-driven bistable overdamped oscillators, governed 
by the Langevin equations 



N 

-y 

4=1 



(1) 



Here £j(t) is a Gaussian white noise with zero mean: 
— 5ij8{t — t')\ e is the coupling constant; 
N is the number of elements in the ensemble, and f(t) is 
a periodic force to be specified later. In the absence of 
periodic force the model (|l|) has been extensively studied 
in the thermodynamic limit N — ► 00. It demonstrates 
an Ising-type phase transition at e = e c from the disor- 
dered state with vanishing mean field X = N Y2i x i 



to the "ferromagnetic" state with a nonzero mean field 
X = ±X §. 

While in the thermodynamic limit the full description 
of the dynamics is possible, for finite system sizes we have 
mainly a qualitative picture: in the ordered phase the 
mean field X switches between the values ±Xq and its 
average vanishes for all couplings. The rate of switchings 
depends on the system size and tends to zero as N — > 00. 

For us, of the main importance is the fact that qualita- 
tively the behavior of the mean field can be represented 
as the noise-induced dynamics in a potential with one 
minimum in the disordered phase (at X = 0) and two 
symmetric minima (at X = ±Xq) in the ordered phase. 
Now applying the ideas of the stochastic resonance, one 
can expect in the bistable case (i.e. in the ordered phase 
for small enough noise or for large enough coupling) a 
resonant-like behavior of the response to a periodic ex- 
ternal force when the intensity of the effective noise is 
changed. Because this intensity is inverse proportional 
to N, we obtain the resonance-like curve of the response 
in dependence of the system size. The main idea be- 
hind the system size resonance is that in finite ensembles 
of noise-driven or chaotic systems the dynamics of the 
mean field can be represented as driven by the effective 
noise whose variance is inverse proportional to the sys- 
tem size H . This idea has been applied to description of 
a transition to collective behavior in jl(| . In [[ll| it was 
demonstrated that the finite-size fluctuations can cause 
a transition that disappears in the thermodynamic limit. 

Before proceeding to a quantitative analytic descrip- 
tion of the phenomenon, we illustrate it with direct nu- 
merical simulations of the model (Q) , with a forcing term 
/(<) = j4cos(Qt). Figure [l] shows the linear response 
function, i.e. the ratio of the spectral component in the 
mean field at frequency Q and the amplitude of forcing A, 
in the limit ^4 — > 0. For a given frequency f2 the depen- 
dence on the system size is a bell-shaped curve, with a 
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pronounced maximum. The dynamics of the mean field 
X(t) is illustrated in Fig. [|, for three different system 
sizes. The resonant dynamics (Fig. ^b) demonstrates a 
typical for stochastic resonance synchrony between the 
driving periodic force and the switchings of the field be- 
tween the two stable positions. 
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FIG. 1: Linear response of the ensemble (|lj)(D = 0.5, e = 2) 
in dependence on the frequency and the system size N. 
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FIG. 2: The time dependence of the mean field in the en- 
semble ^ for D = 0.5, e = 2, A = 0.02, = tt/300, and 
different sizes of the ensemble: (a) N = 80, (b) TV = 35, and 
(c) TV = 15. We also depict the periodic force (its amplitude is 
not in scale) to demonstrate the synchrony of the switchings 
with the forcing in (b). 

To describe the system size resonance analytically, we 
use, following the Gaussian approximation. In this 
approximation one writes Xi = X + Si and assumes 
that Si are independent Gaussian random variables with 
zero mean and the variance M. Assuming furthermore 
that N~ 1 '^ li 6f — M and neglecting the odd moments 
N~ 1 '^ li Si, N~ 1 ^2 i Sf, as well as the correlations be- 
tween Si and Sj, we obtain from ([l]) the equations for 
X and M: 



X 



X-X 3 - MIX + J — V (t) + f(t) , (2) 



M - 3X 2 M - 3M 2 -eM + D 



(3) 



where r\ is the Gaussian white noise having the same 
properties as £i(t). In the thermodynamic limit N — > oo 



the noisy term r) vanishes. If the forcing term is absent 
(/ = 0), the equations coincide with those derived in [0. 
This system of coupled nonlinear equations exhibits a 
pitchfork bifurcation of the equilibrium X = 0, M > 
at e c — 3D. This bifurcation is supercritical for D > 2/3 
in accordance with the exact solution of ([!]) given in [0, 
below we consider only this case. For e > e c the system 
is bistable with two symmetric stable fixed points 



X, 



(2 - e + S)/4 , M = (2 + e - S)/12 , (4) 



(here S = ^f(2 
X = 0, M = (1 - 



e) 2 — 24D) and the unstable point 



£ + v'(l-e) 2 + 12Li)/6. Now, with 
the external noise r\ and with the periodic force f(t) the 
problem reduces to a standard problem in the theory of 
stochastic resonance, i.e. to the problem of the response 
of a noise-driven nonlinear bistable system to an external 
periodic force (because the noise affects only the variable 
X, it does not lead to unphysical negative values of vari- 
ance M, since M is strictly positive at M = 0). 

To obtain an analytical formula, we perform further 
simplification of the system (|2|), ([?]). Near the bifurca- 
tion point we can use the slaving principle to obtain a 
standard noise-driven bistable system 



X = aX - bX 6 




r)(t)+f(t), 



(5) 



where a = 1 + 0.5(e - 1) - 0.5 y/(e - l) 2 + UD, b = 
-0.5 + 1.5(e - l)((e - l) 2 + 12D)- 1 / 2 . A better approx- 
imation valid also beyond a vicinity of the critical point 
can be constructed if we use b = al " 2 instead of 6, where 
the fixed point Xq is given by (^) . Having written the en- 
semble dynamics as a standard noise-driven double-well 
system (||) (cf. 0, ^2|), we can use the analytic formula 
for the linear response R derived in . It reads 
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where s = clNXq/(2D), and T> are the parabolic cylinder 
functions. We compare the theoretical linear response 
function with the numerically obtained one in Fig. 0. 
The qualitative correspondence is good, moreover, the 
maxima of the curves are rather good reproduced with 
the formula (||). 

Above we concentrated on the properties of the linear 
response. Numerical simulations with the finite forcing 
amplitude yielded the results similar to that presented in 
Figs. However, for large amplitudes of forcing (e.g., 
A > 0.1 for £1 = 0.01, D = 0.5, e = 2) a saturation was 
observed: here the response grows monotonically with N. 
This is in full agreement with the corresponding prop- 
erty of the stochastic resonance in double-well systems 
of type ([5]) , where the saturation occurs for small noise 
intensities (cf. Fig. 7 in |ffl]), due to the disappearance of 
multistability for large forcing amplitudes. 



20 40 60 80 

N 

FIG. 3: Comparison of the system-size dependencies of the 
linear response function for frequencies £1 = 0.05 (circles) and 
Q. = 0.1 (squares) with theory (pi). The parameters are D = 1 
and e — e c = 2.5 (where the the exact e c and the approximate 
e c = 3D are used for the ensemble and the Gaussian approx- 
imation, respectively). Inlet: Dependence of the system size 
yielding maximal linear response on the driving frequency SI 
(circles: simulations of the ensemble (Q), line is obtained by 
maximizing the expression (^|)). 

It is instructive to compare the response of the noise- 
driven system with the noise-free case D = 0. With- 
out external force, the ensemble relaxes eventually to a 
steady state solution with some mean field A"; in this 
state each oscillator can be in one of stable steady po- 
sitions of the potential, correspondingly the oscillators 
form one or two clusters. From the clustering it follows 
that the linear response does not depend on the number 
of elements in the ensemble. Our numerical experiments 
demonstrated also that the response is system-size in- 
dependent for large forcing amplitudes as well, where, 
e.g., the force-induced cluster- mergings occur. Thus, the 
effect of system size resonance essentially relies on the 
presence of noise, which breaks the clustering. 

Above we have considered the system of globally cou- 
pled nonlinear oscillators ([!]). The same effect of system 
size resonance can be observed in a lattice with nearest 
neighbors coupling as well. In the thermodynamic limit, 
the Ising-type phase transition occurs in the lattice (if its 
dimension is larger than one). Similar to the globally cou- 
pled ensemble, in finite lattices in the ordered phase the 
switchings between the two stable states of the mean field 
are observed. With the same argumentation as above we 
can conclude that the response of the mean field to a pe- 
riodic forcing can have a maximum at a certain lattice 
size, while all other parameters (noise intensity, coupling 
strength, etc.) are kept constant. We illustrate this in 
Fig. |. 

As the next example we consider the two-dimensional 



FIG. 4: Filled circles: Response of a two-dimensional lattice 
of N with nearest-neighbors coupling for A = 0.02, T — 500, 
D — 0.5, and e = 4. Squares: Response of system (^) (a 
two-dimensional lattice with D = 1.25, e = 30, A — 0.1 and 
T — 140). Circles: the same as squares, but for a globally 
coupled lattice with D = 1, e = 20, A = 0.1 and T = 100. 

nearest neighbor Ising model in the presence of a time- 
dependent external field. The Hamiltonian of the system 
reads 

H = - J SjSj - A cosffif) Si , (7) 

(ij) » 

where J > and s* = ±1. We are interested in the 
dependence of the response of the mean magnetization 
m(t) = jfJ2i s i(t) on tne system size N (for the usual 
stochastic resonance in the Ising model, i.e., for the de- 
pendence of the response on the temperature, see jjj]). 
To calculate the linear response we used the fluctuation- 
dissipation theorem and obtained this quantity by virtue 
of the power spectrum of fluctuations of m(t). The lat- 
ter was found using the Metropolis Monte Carlo method 
on a lattice with helical boundary conditions |l4|]. The 
results presented in Fig. |5| demonstrate the system size 
resonance of the linear response in the two-dimensional 
Ising model. 

As the last example of the system size resonance we 
consider a lattice where each individual element does 
not exhibit bistable noisy dynamics, but such a behav- 
ior appears due to interaction and multiplicative noise. 
This model is described by the set of Langevin equa- 
tions m n 

3 

. (8) 

As has been demonstrated in [ fi"5| , in some region of cou- 
plings e system (||) exhibits the Ising-type transition. If 
an additional additive noise is added to (||) , then one ob- 
serves transitions between these states and the so-called 
double stochastic resonance in the presence of the peri- 
odic forcing |Ti|. As is evident from the considerations 
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FIG. 5: Linear response (in arbitrary units) of the Ising model 
(0) for the temperature T = 2J slightly below the critical 
temperature T c — 2.269J. 

above, such transitions occur even in the absence of the 
additive noise if the system is finite. Thus, the system 
size resonance should be observed in the lattice (|J) as 
well. We confirm this in Fig. [|. 

Another possible field of application of the system 
size resonance is the neuronal dynamics (see, e.g., [[bS[). 
Individual neurons have been demonstrated to exhibit 
stochastic resonance §, While in experiments one 
can easily adjust noise to achieve the maximal sensitivity 
to an external signal, it may be not obvious how this ad- 
justment takes place in nature. The above consideration 
shows, that changing the number of elements in a small 
ensemble of coupled bistable elements to the optimum 
can significantly improve the sensitivity (cf. Q). More- 
over, changing its connectivity and/or coupling strength, 
a neuronal system can tune itself to signals with different 
frequencies. 

Concluding, we have shown that in populations of 
coupled noise-driven elements, exhibiting in the thermo- 
dynamic limit the Ising-type transition, in the ordered 
phase (i.e. for relatively small noise and large coupling) 
the response to a periodic force achieves maximum at a 
certain size of the system. We demonstrated this effect 
for the Ising model, as well as for lattices and globally 
coupled ensembles noisy oscillators. We expect the sys- 
tem size resonance to occur also in purely deterministic 
systems demonstrating the Ising-type transition, e.g. in 
the Miller-Huse coupled map lattice ^0|. The system 
size resonance is described theoretically by reducing the 
dynamics of the mean field to a low-dimensional bistable 
model with an effective noise that is inverse proportional 
to the system size. The stochastic resonance in the mean 
field dynamics then manifests itself as the system size 
resonance. 
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